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l INTRODUCTION 


The free electron laser (FEL) produces coherent radiation 
irem a beam of relativistic free electrons. The principle 
elements of an FEL are the electron beam and a periodic 
magnetic field. The FEL has demonstrated the potential for 
Migh power, high efficiency, and reliability. 

The FEL was first proposed by John Madey in 1970 {1}. 
Since that time, FELsS have been operated and studied by 
institutions throughout the world. Possible uses include 
miGgucumulal,emeGdical, Sciemtific, and military applications. 
The Strategic Defense Initiative Office (SDIO) continues to 
study the FEL extensively for possible use as a land-based, 
ship-based or airborne weapon. 

Chapter II addresses Theater Ballistic Missile Defense 
(PB MD ys. A version of a conventional TBMD is outlined, 
consisting of naval, airborne, and land-based units. The 
future prospects of directed energy, specifically the FEL, are 
also presented. 

Chapter III gives an overview of FEL theory, as background 
for subsequent chapters. Chapter IV describes the Boeing 
Average Power Laser Experiment (APLE), the first FEL project 
with the main goal of producing high average power. APLE is 


in the design stage, and is scheduled to be completed in 1996. 


This research will hopefully provide some additional insight 
into the experiment. 

Chapter V presents the results of research on the proposed 
APLE oscillator design. Numerical simulations are used to 
evaluate both diffraction and pulse effects over many passes 
in the undulator. 

Chapter VI examines both diffraction and pulse effects for 
the APLE amplifier. Here, a numerical simulation that 
evaluated pulse effects over many passes was modified to cover 


only one pass in the undulator. 


ll. THEATER BALLISTIC MISSILE DEFENSE AND THE FEL 


A. BACKGROUND 

Theater Ballistic Missile Defense (TBMD) is one of todays 
Critical national defense issues [2]. The recent Gulf War 
underscored the importance of TBMD, both now and in the 
future. In the current world of only one military superpower, 
the United States can expect its confrontations to be largely 
with countries of the “third world." Many of these countries 
have ballistic missiles, and some are lead by despots willing 
to use them as weapons of terror against the United States or 
her allies. A workable, accurate TBMD system is crucial to 
U.S. interests. This chapter presents a proposed version of 
a system that may work in the near term. This system was 
first conceived and formulated while developing a project for 
a Combat Systems Engineering course (CC4200) at the Naval 
Postgraduate School. Under the guidance of Professor Mike 
Melich, the class studied the concept of TBMD. However, 
G@emectead energy, specifically the FEL, 1s considered the long 
pen SOmurlon tO this ditricult problem. 

Modern ballistic missile warfare has its roots in the 
Second World War. Nazi Germany had limited success using a 
rudimentary ballistic missile, the V-2, to terrorize the 


British citizenry. The British countered ingeniously by using 


deceptive radio news transmissions to draw the missiles off 
Cargew. 

In the recent Gulf War, the Iragis also used ballistic 
missiles for terror. Iraq became aware of the acute political 
impact of employing ballistic missiles during its war with 
Iran [3]. Be waS Go use this lesson to great ettece 1) eee 
firing a total of 72 ballistic missiles [4]. The United 
States and coalition forces countered with the Patriot anti- 
missile system. The Patriot was originally designed for 
protection against manned aircraft, then updated with a 
limited capability against missiles. Politically, the Patriot 
system was a success. It helped to keep Israel out of the 
war, and seemed to protect the fragile coalition itself. 

The military success of Patriot is now in question, 
however. As for the Patriot’s overall track record, and for 
its value in theater missile defense, there seems to be little 
agreement. Nevertheless, experts agree that TBMD is extremely 
important, spurring much recent research and development. 
According to a Strategic Defense Initiative Organization 
(SDIO) official, funding for this R&D 1s expected teubevaea 
$500 million per year through fiscal year 1993, and to peak at 
about $1.51) 1] tense eS ela 

In the near term, the most viable candidates for active 
defense are the "missile-on-missile" engagement systems. 
However, according to Major S.A. Fleet, USAF, of the Air 


Command and Staff College, “one of the most difficult 


challenges facing weapons designers today is how to build a 
missile that is agile and accurate enough to destroy attacking 


warheads" [6]. 


B. A NEAR TERM TBMD 

The Navy can be a major player in Theater Ballistic 
Missile Defense. A near term TBMD system would be built 
around the Aegis cruiser, the Navy’s most capable anti-air 
warfare (AAW) ship. This system would provide a mobile, quick 
response TBMD for the protection of U.S. interests ashore and 
at sea. This recommended defense system 1s purely 
hypothetical; any association with any other conceptual 
systems which may or may not exist is totally coincidental. 

The heart of this TBMD system is an updated Aegis cruiser. 
This “update" would be a variant of the most capable Aegis 
phased array radar, the SPY-1D. The update would include 
specific hardware and software changes to tailor Aegis to the 
TBMD role. Battle groups will deploy with at least one of 
these ships. The Aegis cruiser performing this role can be 
called "SLINGSHOT", an acronym for Surface Launching Ship on 
Die aiee re 

Defense Support Program (DSP) satellites provide the first 
ingacabion Sof, a ballistic missile launch. .[5)}. These 
satellites were originally designed to detect the relatively 
large plumes of intercontinental ballistic missiles. 


Detection algorithms and processing software in the satellites 


ground stations were updated during the Gulf War, providing 
good capability against the relatively short-range theater 
ballistic missiles [5]. These reconnaissance satellites, with 
their infrared sensors, provide an initial warning and a very 
rough estimate of the intended target. The ground-based 
processing stations relay the warning to SLINGSHOT via 
satellite communication links. 

The Aegis missile battery is the first level of defense 
for coastal areas. Upon receipt of a cue fromeythemeee 
satellite, SLINGSHOT will search for the missile using its 
Special TBMD radar mode (range on the order of hundreds of 
miles, altitude on the order of several hundred thousand 
feet). After detecting the ballistic missile and establishing 
a firm radar “track," SLINGSHOT will engage the target with 
long range TBMD missiles. The Aegis system will then perform 
kill assessment, the results will determine whether the 
incoming target should be re-engaged. 

For areas further inland, fighters from the battle group’s 
aircraft carrier will be the first level of defense. The 
cycle of the aircraft will be like a "stack." Two pairs of 
F/A-18 will form an upper and lower racetrack at 35 kft and 25 
kft in the vicinity of the high value area (HVA). The 
altitudes will largely depend on weather and aircraft loading. 
Meanwhile, one KA-6 will loiter in the HVA to provide fuel for 
the fighters. There will be another tanker (S-3 or F/A-18 


with buddy stores) available between the HVA and the carrier 


CO act as back-up. «On deck, two F/A-18 and one tanker are at 
alert 5 (five minute warning prior to take-off) and two F/A-18 
tieeeat alert 152 When the upper CAP’s cycle is complete, they 
return to the carrier (“homeplate"). The lower CAP 
Simultaneously moves to assume the upper racetrack position. 
The alert 5 aircraft launch and assume the lower CAP position. 
This “stack” continues in this manner, with the fighters low 
on the stack moving progressively up until they are the upper 
CAP. Variations of this include the CAP doing double cycles 
with tanking in between. This would decrease the burden on 
the carrier. 

A broad view of a possible TBMD air engagement will 
illustrate this tactic. Assuming a 600 km ballistic missile 
flight path, one can expect a terminal speed of about 5-7 mach 
and a terminal angle of about 48-50 degrees (from horizontal). 
Positioning the CAP on the far side of the HVA provides 
adequate horizontal space for maneuvering into’ launch 
parameters. Upon receiving a cue from the DSP satellite, 
SLINGSHOT will relay the cue to the upper/lower CAP via a high 
Speed umaacta slink. The upper CAP (loitering subsonic to 
conserve fuel) now turns to the target bearing, and maneuvers 
to meet launch parameters. These maneuvers should take 
several minutes. The pilot closes his firing key when he 
meets all parameters; SLINGSHOT can now fire the missile 
remotely via data link. The other upper CAP aircraft follows 


in trail and conducts a similar engagement, thereby increasing 


thes. pmebabimlity sof) ki) eae. Immediately after the 
engagement, these CAP will tank with the KA-6 and then proceed 
to homeplate. Meanwhile, the lower CAP have moved into the 
upper position, and the former alert 5 aircraft are now in the 
lower CAP position. Aircraft control will occur via several 
voice communication circuits and the aforementioned data link. 
The carrier uses launch and recovery to handle the aircraft 
during launch and recovery operations. SLINGSHOT controls the 
CAP on the war net, while tanking is controlled by the carrier 
on the admin net. [7] 

The Patriot system provides the innermost layer of 
ballistic missile defense. However, since Patriot is 
relatively immobile, it may take an additional five to ten 
days to be on station and operational. The Patriot battery 
may be "linked" with the other defense systems or it may stand 
alone. 

The Navy does not have any anti-ballistic missiles in its 
current inventory. SM-2 block 4 and Erint (Extended Range 
Intergepror ) are among the candidates for shipboard 
deployment. SM-2 block 4 is a future update to the Standard 
missile family; SM-2 block 2 is the latest available variant. 
Erint, a product of LTV Corporation, is designed speciitieaiee 
for TBMD. It is a small, agile "“hit-to-kill" missile with 180 
small maneuvering rocket motors [5] that supplement its fins: 
Erint is a good candidate for the missile to be deployed on 


aircraft due to its relatively small size. 


C. THE MOTIVATION FOR DIRECTED ENERGY WEAPONS 

The advent of gunpowder, followed by Napoleonic warfare, 
and Nuclear warfare were the great "watersheds" in the history 
of warfare [8]. It is proposed here that the use of directed 
energy weapons will be the next great "watershed." This 
natural evolution will occur because "missile-on-missile" 
systems have nearly reached the end of their utility. 

In this section, the author makes several assumptions. 
The range of a surface-based directed energy weapon is assumed 
to be on the order of tens of kilometers; an airborne laser is 
assumed to have a range of hundreds of kilometers. Weather 
effects, which could be significant, are not considered. 

The speed of modern missiles places a high premium on 
reaction time. An incoming Theater Ballistic Missile (TBM) 
with a speed of 5-7 mach (1.7-2.4 km/s) and a range of 600 km 
leaves about 6-8 minutes for detection and engagement. as 
kill assessment indicates a miss, there may be no time left 
for a re-engagement. The high speed of TBMs also complicate 
the £Luzing problem in proximity warheads, requiring extremely 
BaioEwbaZing aia Dlast action. High speed targets also place 
additional constraints on hit-to-kill anti-ballistic missiles 
(ABM), requiring them to be both extremely fast and agile. 
The conventional ABM fly-out time must also be considered. A 
mach 4 ABM requires ~40 seconds to travel to a 55 km (30 NM) 
intercept, while a directed energy weapon requires only 


0.00018 seconds. However, for a laser, dwell time on target 


must be consideredyenequiningsacci tional time wonmeenerc: demson 
a Second. 

Magazine depth, the amount of available ordnance, could be 
another critical issue. If the enemy fires a large number of 
TBMs during a short period of time, battle group assets could 
rapidly become depleted. However, some directed energy 
systems could conceivably have an "unlimited magazine." 

In the first or second decade of the 21** century, 
directed energy weapons could replace the conventional TBMD 
weapons described in this chapter’s previous section. The 
Shipboard high energy laser (HEL) could replace the surface- 
to-air missile (SAM). The expected range of this HEL would be 
on the order of tens of kilometers. An airborne HEL, based on 
a relatively large airframe, would replace the TBMD fighters. 
The airborne laser’s range would be on the order of several 
hundred kilometers. 

Range is dependent on the level of radiation intensity on 
target mrequincadmror a Nie. This level, termed target 
hardness, differs with each type of missile target. According 
to Major Gary Danczyk, USA, of Los Alamos National Laboratory 
(LANL), a nominal TBM would require »>1.0 kJ/com? (absorbed) 
over a@ minimum axial length to ensure unstable crack 
propagation. Using a dwell time of 1 second, an absorbed 
intensity of >1.0 kW/cm* is required. 

Range is also critically dependent on atmospheric effects. 


The total beam attenuation depends on the atmospheric 
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conditions, wavelength of the radiation, and the total 
atmospheric path length. Losses, particularly atmospheric 
Scattering and absorption may be significant, raising the 
power requirements of the laser. Absorption can lead to an 
effect called "thermal blooming," a deflection and defocusing 
of the beam. These atmospheric effects can be minimized by 
using laser wavelengths that give a relatively large percent 


transmission. Here, a tunable laser could be advantageous. 


[9] 


D. WHY THE FEL? 

The FEL is the directed energy weapon of choice for the 
Qe weecentury., It is all electric, with a potential for high 
efficiency, an "infinite magazine,“ and a continuously tunable 
wavelength. 

VADM Kihune, USN, Vice Chief of Naval Operations for 
Surface Warfare, ina brief at the Naval Postgraduate School, 
addressed the goal of having the "all electric ship" by the 
year ~2015. This goal is linked directly to supporting the 
shipboard HEL. The FEL and the all electric ship are an 
excellent match because of the FEL’s requirement for large 
electric power. 

Another advantage of the FEL is its capability for 
relatively high efficiency. The wallplug efficiency of an FEL 
is the ratio of the average radiation power output to the 


electric power input. The single-pass extraction efficiency, 


ital 


Nn, is the fraction of the electron beam power that is 
converted to optical power over a single pass. The extraction 
efficiency of the APLE amplifier is estimated to be ~5%. An 
extraction efficiency of n~40% has been demonstrated. Some 
FELs, such as the CEBAF experiment, estimate wa awal lolug 
efficiency as high as 40% {10}. 

Magazine depth is an important issue, as mentioned 
previously. The FEL “magazine" is likely “deeper" than that 
of other lasers. This is because the basis of the FEL is only 
electrons and a magnetic field. Other lasers use some type of 
active medium which would likely require pervedte filcering 
and replenishment. 

A major advantage of the FEL is its readily adjustable 
wavelength. Conventional lasers do not have this capability 
because the output wavelength is a function Of €hemnaminea 
resonant frequency of the atom or molecule in the active 
medium. The FEL wavelength, A, can be changed by changing 
the electron beam energy or the undulator field strength (see 
equation 3-8). FELs have demonstrated operation from 240 nm 
to 9 mm, and have exhibited continuous tunability of a single 
FEL over the operating wavelength by a factor of 6 {11]}). This 
tunability would enable the FEL to take advantage of the 


prevailing atmospheric conditions during any TBMD engagement. 


2 


lil. FREE ELECTRON LASER THEORY 


A. BASIC FREE ELECTRON LASER PHYSICS 

The FEL consists of two basic parts, an accelerator and an 
Gneulators or  “wiggler a: The accelerator produces a 
relativistic electron beam, while the undulator produces a 
periodic magnetic field ( see Figure 3-1) [12]. The electron 
beam oscillates in the transverse direction as it moves along 
the undulator. The periodic, transverse acceleration of the 
electrons causes them to radiate ina forward cone. Some of 
this spontaneous radiation can be stored ina laser resonator, 
formed by placing two concave mirrors (one semi-transparent) 
at either end of the undulator. The coupling of the optical 
mode and the electron beam in the undulator enables kinetic 
energy to be extracted from the electrons. This kinetic 
energy is converted to electromagnetic radiation, leading to 


stimulated emission and coherent radiation. The FEL mechanism 


electron beam 


undulator 


¢ 4 
wre wos > 
ae 





resonator mirror 


Figure 3-l. The FEL schematic. 
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is fundamentally classical and may be described with 
electromagqnetiwesrncon, | 138 
The accelerator provides the relativistic electron beam 


for use in the undulator. The beam energy, ymc*, ranges from 


a few MeV to a few GeV. The Lorentz factor, y=1//1-8?, where 


§=/c is typically in the range of y=10* to y=10°. The 
typical FEL uses a beam energy of about 50 MeV with y=10°. 
The current ranges from 1 A to 1000 A, but most FELs fall into 
the range from I=10 A to 100 A. The typical power carried by 
the electron beam is therefore =5 GW, so that even for small 
efficiency the FEL can yield high peak power. However, since 
the electron beam from an accelerator is usually pulsed, the 
average power is usually much smaller. 

The undulator consists of a periodic magnetic field that 
1s either circularly or linearly polarized. The undulator 
length, L, can range from 1 m to mone than 25 mm) Dive 
typically 5 min length. Each undulator period, 4, is abeus 
5 cm long, but can range from 1 cm tor l0"cm. The resulemig 
number of periods, N, 15 typically  N=l00, 3 The streng emia 
the undulator can be described by the undulator parameter, 
K=eBA,/2nmc*, where e=|e| is the charge magnitude of an 
electron, mis the mass of an electron, and c is the speed of 


light. The rms magnetic tield strenqmiye 2/2 15 typwesme, 


a few kilogauss so that K=1. 


14 


Figure 3-2 illustrates the interaction of an electron beam 
with the electromagnetic radiation within an FEL resonator. 
At the top frame, electrons enter the undulator where they 
encounter the periodic magnetic field, and begin to oscillate 
in the transverse direction. This oscillation causes the 
Siecrenons tO emit light (photons) in a forward cone. The 
middle frame shows a single period within the undulator. The 
radiation moves through one period of the undulator magnetic 
field as the slower electron follows a sinusoidal path. The 
bottom frame shows that the effect of the optical field on the 
electron depends on the relative phase between the electron 
and the optical field. Here, the electron sees the maximum 
radiation electric field, oriented in the same direction as 
the transverse motion of the electron. Since the electron is 
negatively charged, this field exerts a retarding force on it. 
The electron consequently decelerates, giving up energy to the 
radiation field. The energy given up by the electron is 
transferred to the laser beam, thereby amplifying the 
Paci at lon. Energy may also be absorbed by the electron, 
causing it to accelerate and resulting in a loss in the 
optical field. Therefore, in a beam of electrons extending 
over many wavelengths of light, some gain energy and some lose 


energy. This results in "bunching” of the electron beam [14]. 
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FEL optical field-electron beam interaction. 
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B. ELECTRON DYNAMICS AND THE PENDULUM EQUATION 
The fundamental feature of the FEL is the axial bunching 
described in the previous section. Assuming a plane wave, the 


optical fields are [15] 


E,=E(cosy,-siny,0),  8,=E(siny,cosy,0), seal 


s 
where E£— is the optical field strength, WkzZ-wt+¢@ with carrier 
frequency wekc, k=2n/A is the optical wavenumber correlated 
with the wavelength A, and ¢ is the optical phase. 

The electrons periodically deflect in the transverse 
direction as they travel along the undulator axis. The ideal 


form of the helical undulator field is 


B=Blcos(k,z),sin(k,z),0] , 3-2 


where B is the peak undulator field strength, k,=2n/A is the 
undulator wavenumber associated with the wavelength 4, and z 
Peomedewairection along the undulator axis. 

The interaction of the electrons and the optical fields 


inside the undulator is described by the relativistic Lorentz 


tomec pls] 
—S (yB) =-—& (H+ (Bx8)1 , 3-3 
dt mc 
and 
LCN oerequepl ag fais : 
dt ale aa 


Using the combined fields of the undulator and radiation, the 


transverse motion with perfect injection can be found to be 
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B, = -F [cos (752) 7 S11) kee - [siny,cosy,0], 3-5 


where A=eEA/2nmc* is related to the optical vector potential 
and is analogous to the definition of K. 
Using equations 3-3 and 3-4, the change of electron energy 


is found to be [16] 
2 Che 
=e OS ( (+d) =, 356 
“ y a Ss (C+) 


where ¢€ is a dimensionless variable that describes the 
electron’s phase with respect to the combined optical and 
undulator fields. Electrons with phases such tha, 
-m/2<(C+¢)<nm/2, will absorb energy from the optical fae 
CauSing them to accelerate. The electrons with phases, 
m/2<(€+¢)<3n/2, will lose energy to the optical field causing 
amplification. Given a random distribution of a large number 
of initial phases, about half the electrons gain energy and 
about half lose energy. However, as described earlier, under 
the influence of combined fields, the electrons can be made to 
"bunch" at desirable phases, causing gain during a pass down 
Phe Une ul aieeme 

To further clarity thes dynamtesmeet electrons in ee 


undulator, it is convenient to rewrite 3-6 in terms of ¢ alone 


for low efficiency. Using the relation y~?=1-f*-f%, and the 


conditions y>>1l, 8,=1 and ®*=K2/y2, it can be shown that the 
y, z 1 7 
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individual electrons follow motion described by the simple 


pendulum equation [16] 


t =v=|lalcos (¢+6) ald 


where v=(=L[(k+k,)B,-K] is the dimensionless electron phase 


velocity, and |a|=4nNeEKL/yjmc* is the dimensionless optical 


mrela Strength. In 3-7, "°" means d/dt where t=ct/L is the 
dimensionless time; a pass through the undulator is described 
by t=07l1. The electron phase velocity, v, 1S a measure of the 
resonance between the electron beam, the undulator, and the 
optical’ fmeld. 

When v=0, the FEL is said to be at "resonance." The 


resonant wavelength is 


This equation shows that the operating wavelength of an FEL 
Game changed by changing the undulator wavelength 4, the 
ie uau@@um@emmaneter K,wor the sewectron beam energy ymc*. 

For high efficiency, when the change in phase velocity, v, 


is comparable to N, the pendulum equation can be written more 


ae@eunarely as 


Y 
1| 
<o 


=|a|(1-3v/4nN) cos (C+o) . At 


i 


C. THE FEL OPTICAL WAVE EQUATION 
In developing the optical wave equation it 1S appropriate 
to begin with the Maxwell’s wave equation in the Coulomb gauge 


[15], 








where, W=+0+0, A=(E/k) {[sinp,cosy,0] is the vector 


potential representing radiation trom a helical undularer mara 


— 


J, is the? transverse Current density ~paome Chemeranc emo 


motion of the electron beam. The optical field envelope is 


taken to vary slowly in time during an optical period (EX <a, 


d<<wh), and in space over an optical wavelength 
(EC CLE Ce Kida a. This is the "slow-varying amplitude and 
phase approximation." Using the above approximation, the left 


side of 3-10 can be written as 


Aaa Sees 4 [cosy, -siny, 0] + 





Me ee oc 
v C2 6z cot 





-25| 58 += 4 [siny,cosy, 0] 


Q) 
Ct 


Replacing J, with the sum of all single-particle currents, and 


assuming low efficiency so that the changes in y are small 
even at saturation, and averaging over an optical wavelencen, 


the wave equation may be expressed as 
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ae tS ae|zee= - 28 oii ae 3-12 
MHere S@—@:,tk)\z—-Wt, po is the electron particle density, ,and 
<...> represents the ensemble average over all the electrons. 
In the special case where the electron pulse is long, it has 


no spatial dependence. The wave equation is then independent 


of z and becomes 


eek ee 3-13 
where a=|ale’® is the complex dimensionless field, and j is 


the dimensionless current density 


a= ARNeKLEe** fe 8N(enKL)2p 
y?mc? ; y?mc? 


3-14 


For a linear undulator the derivation is similar, except 
that K is replaced by K(J)(§)-4,(&)), where E=K=/2(1+K*). The 
Bessel function factor accounts for reduced coupling in the 
Vinear wundulator due to fast periodic Zz motion of the 
Pre cwmoisa ta) | . 

For the high efficiency case when the changes in v can be 
Bewetacrge was N, the wave equation can be more accurately 


written as 


8=J<(1-v/4nNn) et . 3-15 


Ok 


D. FEL DIMENSIONLESS PARAMETERS 

A few dimensionless variables can summarize recurring 
combinations of physical parameters in a complex problem like 
Ice) aby several dimensionless variables are used to gain 
insight into relevant physical processes without detailed 
calculations or simulatioms pi. 

Some dimensionless variables have already been introduced. 
The undulator parameter, K, 1S an important quantity 
characterizing many FEL properties such as the electron 
deflection angles. The electron phase and phase velocity, ¢ 
and v, define bunching of electrons and their bunching rate. 
The optical field strength |a|[, can also be used as an 


indicator of expected electron bunching rate since y-«la]. 


The dimensionless current density j (see 3-14) determines 
the response of the optical field to electron beam bunching. 
It also provides the coupling between the electron beam and 
the light wave. The FEL’s gain regime 1s dererminedes, sene 
value of j: when j-l, or less, the FEL gain is low, and when 
j>>1 the FEL gain.is high. 

The slippage distance, NA , is used to normalize 
longitudinal distances in the FEL. The slippage distance is 
the distance that light passes over an electron as the 
electron travels = through ) encmeuneewaven: It as” tive 
characteristic distance in the trapped-particle instability 


and short-pulse evolution. At resonance, exactly —one 


ae 


wavelength of light passes over an electron as the electron 
travels through an undulator wavelength. 

iicmcharacterist temmmapta cal ibeam radius, athA;m i, 
normalizes transverse dimensions in the FEL. This is done by 
multiplying the dimension, such as the electron beam radius, 
Bymanmnei’ . 

The "filling factor", F, is the ratio of the electron beam 
area to the optical mode area, averaged over the undulator 
Tengu, and sadetermimess the effect of transverse optical mode 


CoiSteextacons 


E. THE HIGH GAIN REGIME 

When the dimensionless current density is large enough so 
that j>>1, the optical wave amplitude and phase change 
Significantly during the FEL interaction, usually resulting in 
high gain. Maximizing j with the aim of achieving high gain 
depends on several factors. At short wavelengths, j tends to 
be relatively small because more relativistic electrons are 
difficult to bunch (refer to equations 3-8 and 3-14). Working 
at longer wavelengths will increase j, but the optical 
wavelength is usually determined by the laser’s application. 
The length of the undulator, L, can also be increased, but 
this increases the FEL’s sensitivity to beam quality. On the 
other hand, when the current increases markedly, the beam 
quality from the accelerator starts to decrease. In other 


words, the penalty for achieving a large j is often a 
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corresponding decrease in beam quality leading to poor gain. 
However, this degradation can be offset by mode distortion 
(see section III.I.). 

The wave and pendulum equations are valid in the high gain 
regime, and can be solved numerically to explore phase space 
evolution. In weak optical fields, the equations can also be 
solved analytically. Expanding the pendulum and wave 
equations, and integrating over all initial phases ranger we 
can obtain an integro-differential equation governing the 


evolution of@ehe opticalgineld | ie], 


ce) 


ala =a) [ de’t! F(t’) SEE Rt, 3-16 
2 Jo 
where F(t') =[dgf(q)e4* is the characteristic function of the 


distribution “f(¢q), f(q) is the distribution of eniteua 


electron phase velocities vV.—V7 cuaootr same [ea £1 ¢) 


The integral equation 3-16 is valid for both high and low 
current. The equation can be solved analytically or 
numerically. In. the limit of weak -optical fields, the 
integral equation has an advantage over a simulation that must 
follow the evolution of many sample electrons. In the high 


current limit j>>1, the integral equation can be _ solved 
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analytically to obtain the gain for a perfect beam [12], 
G(t) = te (F/2Ovie 3-17 
2 


At™high current, intér—-electron Coulomb fomees may effect 
SmccUron DUNCHING ——slUlesce Perececuacemnegligible except, for a 
value of 7 that would result in extremely high gain; in this 
case, the high gain itself would already cause much more 
Gmandeve COmeecrivG =erreecus."secoumomb forces can also cause 
distortion of the electron pulse, but this is usually an 


tMmomgniiicance effect (12). 


F. THE TAPERED UNDULATOR AND HIGH EFFICIENCY 

High efficiency is usually desirable in the FEL. A 
tapered undulator is required in order to achieve high 
CS erdeulommerrticlency, 7, the traction of electron beam energy 
converted to light. The maximum extraction efficiency of an 
untapered undulator is n=l1/N and is typically only a few 
percent [19]. In a tapered undulator, about 50% of the 
electrons are coherently decelerated in ponderomotive 
potential wells, or “buckets," formed by the combined 
undulator and radiation fields. These "trapped" electrons 
continue to lose energy to the light wave as the undulator 
length is increased, so that nm may exceed =1/N. 

The beginning section of the undulator should be untapered 
with a constant peak magnetic field B and spatial period A. 


This allows the power in the light wave to grow to saturation 
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where trapping begins. After saturation without taper, the 
laser power will decrease with undulator length as electrons 
absorb power from the light. Equation 3-8 shows that as the 
electron energy 7 decreasesvatisaturation,weither A woneiaean 
be changed to compensate, maintaining the resonant wavelength 
through the undulator. The basic.ddeanof tapering susmiee 
change sthe Gundulaser,) pasameters, B on Aco Sitcoms 
resonance with the radiation field as the trapped electrons 
lose energy. 


Decreasing the undulator field strength, or the undulator 


wavelength, results in a phase acceleration, C= 8+... where 


6=87INN ges: red ANG Wane nar is the desired efficiency assuming 
50% of the electrons are trapped. When N>>1, the phase 
acceleration can be significant with only a small change in 
the undulator properties. 

Recall that the pendulum and wave equations, 3-9 and 3-15, 


had a correction factor for  iigqueee: 2 vem. Forma high 


efficiency, tapered undulator y=§, and av=6at=6 for at=l. 


Using the identity in the previous paragraph with a desired 
efficiency of 5%, yields 6=4nN/10. Taking <av>=4nN/10 and the 
fact that V>AV, results sti) COrEect 20m factors ope 
1-3<v>/4nN=0.7 and 1-<v>/4nN=0.9 in 3-9 and ci lie 
respectively. On the other hand, av=n for the untapered 
undulator, resulting in the correction being =1 and therefore 


not used. 
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G. SHORT OPTICAL PULSES 

Short, picosecond long electron pulses from an RF 
accelerator drive many FEL oscillators and some FEL 
amplifiers. When the electron pulse length is comparable to 
the slippage distance NA, short-pulses will effect the FEL 
interaction [20-25]. feenemweleceron pulse length is much 
longer than the slippage distance, the pulse is considered 
onc The performance of the APLE oscillator is somewhat 
independent of the short optical pulse effects, but the pulse 
effects that do occur will effect the amplifier. 

Within the FEL oscillator short electron pulses lead to 
short optical pulses via spontaneous emission. Thi smeoht 
bounces inside the resonator between concave mirrors separated 
by a distance S. The electron pulses from the accelerator 
must be synchronized to arrive at the beginning of the 
Pndimlaton c—0; coincident ewith the rebounding light. The 
distance between the electron and the optical pulse at t=0 is 
termed "desynchronism," d=-aS/NA , and is normalized to the 
slippage distance. When d=0, there is exact synchronism 
between the light and the electron pulse. 

Recall that on each pass down the undulator, the unbunched 
electron pulse starts at t=0 and slips back with respect to 
the light as bunching develops. The bunched electron pulse 
tends to drive the trailing edge of the light pulse, 
distorting it on each pass. As a result, the centroid of the 


light pulse slows to below c. To compensate for the slower 


2 


speed of the light pulse, the path S must be reduced sovthat 
d=0.00120.1 [12]. In practice, moving one of the resonator 
mirrors a small distance aS adjusts the value of d. The 
adjustment must be chosen carefully since the incorrect value 
of d will adversely effect the overlap of the electron beam 
and the optical mode, thereby reducing gain. 

The general features of short-pulse behavior have been 
observed in many simulations and several experiments [12]. At 
small desynchronism, d-~0 or greater, the power is often large 
enough to cause the trapped-particle instability, a broad 
optical spectrum, and a broad electron spectrum. The optical 
pulse is short and is centered on the electron pulse. The FEL 
tends to be unstable at small values of desynchronism, since 
a small change in d can make a large difference in steady- 
state power. 

For “Farger cd, thes FEl becones se nionmemorante. At large 
desynchronism, the steady-state power is smaller due to the 
reduced coupling, and the trapped-particle instability usually 
does™not occur™’ The fimal optical pWise cangbe much lenge. 
than the electron pulse, resulting in a narrow power spectrum. 
The electron spectrum is narrow in the weaker optical fields, 
and the center of the optical pulse may actually be well ahead 
of the electron pulse. The FEL’sS operating range in d is 
typically ad=0.1, so that the resonator length must be 
adjusted within a surprisingly small range aS=-5 um for N=10° 


and A=lum. 
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H. THE FEL TRAPPED-PARTICLE INSTABILITY 

When either the FEL amplifier or oscillator reaches power 
levels beyond the onset of saturation, the FEL may experience 
wiew cLapped—-particlé instability. Electrons can become 
trapped in deep potential wells in phase space when under the 
influence of strong optical fields |a|>>m. The oscillation of 
the beam current that is trapped in these wells can drive the 
carrier wave unstable, causing sideband frequencies to grow 
from noise. The trapped-particle instability can then modify 
or destroy the coherence that is initially established in weak 
optical fields. This instability can be exhibited in a wide 
mane c OUsUNdUt ater Aesignus operating in strong fields [26-31]. 

The general features of the trapped-particle instability 
fimeicemmemascriiator depend only On the dimensionless current 
j and the loss factor Q. For sufficiently low values of j and 
QO, the trapped-particle instability does not occur; however, 
ineuweasing either j er QO can result in the instability. 

The trapped-particle instability in an FEL is a mixed 
blessing. A high-power laser with a narrow spectrum and no 
sidebands is often desirable to the experimenter. However, 
miemmemesence. of sidebands indicates the achievement of high 
POWemeand@enat Ghe FEL 1s G@pemating well.» Inefact, the FEL 
attains even higher power and efficiency as more sidebands 
develop. High power tends to result in both a broad electron 
spectrum and a broad optical spectrum, but tapering tends to 


diminish the instability. 


Zo 


In the amplifier, the eleéetron beam intenacts (wien tne 
optical field for only one pass. Therefore, only the largest 
values of j=10° have a possibility of developing a significant 
Sideband. Because of this requirement, the trapped-particle 


instability has never been observed in the FEL amplifier. 


I. FEL MODE DISTORTION 

The FEL interaction modifies both the optical phase and 
the wavefront. When considering transverse effects, the 
optical field is given transverse spatial dependence, 


a(x,y,t). In this case, the parabolic wave equation is 


3-18 


gs Te) aae 


e Vite a (ogy 1) = yeu ae, 


where Wi=0,+0 must be used “co Winclude” themsetfeccem—a: 


diffraction [12]. In either strong or weak fields, theweas 
interaction causes a phase shift that is opposite to that of 
natural diffraction, and therefore focusses the light back 
into the electron beam along z. For large current j and a 
small electron beam, this focussing can decrease the 
transverse area of the optical mode, thereby increasing the 
effective F. The net result is an increase in the effective 


dimensionless current, jE, Yeading, toubrgher gain. 
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J. BEAM QUALITY 

Recall Chat eptimizingean FEL experiment often leads to a 
design trade-off between high beam current and poor beam 
quality as defined by the accelerator. A given application 
usually fixes the optical wavelength A so that the 
dimensionless current, j«xIN°A’’*, must be maximized with the 
current I and the number of undulator periods N. However, 
increasing I tends to decrease beam quality from the 
accelerator, while increasing N increases the FEL’s 
sensitivity to beam quality because the gain bandwidth is 
narrower. The consequence of poor beam quality is a reduction 
in FEL performance due to degraded bunching [12]. 

Beam quality is generally described by two qualities: the 
beam energy spread, ay/y, and the beam emittance, e. The 
energy spread is a measure of the spread in electron velocity 
magnitudes. The Gaussian distribution of energies with rms 
emcdorAy/y, can be characterized by o,=4NNay/y; when o,=r 
Dimeiingu:s impaired. The value of o. gives the spread in 
electron phase velocities av in the beam. Emittance is a 
measure of the spread in electron velocity directions or 
angles. 

Emittance is given by e=r§8 where fF is the average beam 
radius, and § is the average beam angular divergence. When 
the beam is focused to a small radius f, the angular spread § 
increases keeping the emittance e« fixed. As the beam is 


2 


accelerated, and the beam energy ymc is increased 
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Significantly, the normalized emittance e .=ye=yr8 tends to 
remain constant and is a useful measure of beam quality. A 
Gaussian spread in angles with rms value 6, results in an 
exponential Grstribueven of phase velocities with 
characteristic width og=4nNy°O0°/ (1+K*). When og is equal to 


or larger than nm, electron bunching is impaired. 
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IV. THE AVERAGE POWER LASER EXPERIMENT 


The Average Power Laser Experiment (APLE) is one of the 
most significant studies in the free electron laser field. 
Although the FEL has potential for use in a wide range of 
disciplines, it currently has been designed for only mediocre 
average power. This experiment, now under development, 
intends to demonstrate the high-average power capability for 
the first time. 

APLE [32] 1S a cooperative effort between Boeing 
Corporation and the Los Alamos National Laboratory (LANL). 
Its objective is to design, construct, and demonstrate an 
integrated, hagh duty factor, 10 micron wavelength free 
electron laser by 1996. The primary design goal is to sustain 
100 kW average power for at least three minutes. The design 
parameters are listed in Table 4-1 [33]. The device will 
operate in the single accelerator, master-oscillator, power 
amplifier (SAMOPA) mode, with the oscillator providing a seed 
laser of about 100 watts to drive the amplifier. In the 
acronym SAMOPA, "SA" refers to the APLE accelerator that 
powers botn the oscillator and amplifier; “MO“ refers to the 
FEL oscillator, and "PA" refers to the FEL amplifier. The 


cost of APLE is approximately 100 million dollars; initially 
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TABLE 4-1. APLE FEL DESIGN PARAMETERS. 


Configuration - Single accelerator oscillator amplifier 
Laser power - 100 kW 
Wavelength - 10.6 pm 


Run time - 3 minutes 


Electron beam Oscillator Amplifier 


17_Mev 34 MeV 


AVG] CUEBrent Oe23 # 0-23 & 
(macro) 


Avg. power (25% 1 MW 2 MW 

duty) 

Pulse length 10 ms 10 ms 

(macro) 

Pulse rate 27 MHz Zo Miz 
(Eee) 

Pulse length 18 ps 
peer oem. 

peak current si current oo A i A 
Normalized 80m (mm-mrad) 80m (mm-mrad) 
emittance 

Energy spread 0.3% 0.75% 

(micro) 
a 
Extraction 0.01% 

| ie, 


Laser [Laser power si 10 ol W 100 ho Gael 
[Resonator | 22 m concentric [== 
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phieweuad date was seumar 1994) but the project will run longer 
daeeeo runding limitations 1/34]: 

The existing Boeing FEL facility in Seattle, Washington, 
is the sight of the experiment. iges laboratory complex 
provides virtually all the required housekeeping support. A 
visit to Boeing in November 1991 revealed a large two story 
building with a basically empty second floor; this area will 
eventually house the laser. The amplifier wiggler and many 
crates of support equipment are on location. Much of the 
ieaawane and SOLtware to be usedvon APLE are on site because 
of use on previous FEL experiments. 

The main components of APLE are shown in Figure 4-1 [32]. 
They are the photoelectric injector, the accelerator, the 
oscillator and amplifier wigglers, the beam transport optics, 
and the electron beam dump. The accelerator is separated into 
two sections by the oscillator wiggler and associated 
hardware. Each accelerator section 1S powered by two 
klystrons. 

A general system description starts with the production of 
electron micropulses in thle photoelectric injector. The 
electrons exiting the injector have a normalized emittance 
(90% edge) of 80n mm-mrad. These electrons are then 
accelerated to 17 MeV at a current of 0.23 A by the first 
section of the accelerator. At the end of the first 
accelerator section, the beam is focused, bent, then re- 


focused prior to entrance into the oscillator wiggler. The 
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The APLE layout. 


Figure 4-1. 
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electron beam exiting the oscillator wiggler is then bent into 
the second section of the accelerator. The energy spread (90% 
edge) in the wiggler is 0.3%. Two interleaved chicane 
bunchers separate this section into two segments. The first 
segment boosts the electron energy to 25.5 MeV and is phased 
to impart a 2.6% energy slew on the electron pulses. This 
slew 1s necessary in order to bunch the beam in the subsequent 
chicane bunchers, which compress the beam micropulses from 60 
ps to 18 ps FWHM. The reduction of the electron micropulse 
length from 1.8 cm to 0.54 cm increases the peak current in 
order to increase the amplifier gain. The final accelerator 
segment increases the energy to 34 MeV at a current of 0.23 A 
and removes the energy slew from the beam before it enters the 
amplifier wiggler. The energy spread in the wiggler is 0.75%. 
After exiting the amplifier wiggler, a magnet directs the beam 
into the main beam dump. The built-in flexibility of the APLE 
layout allows the electron beam to travel straight-ahead, 
through the oscillator wiggler only, or through both the 
oscillator and amplifier wigglers. Preconditioning of the 
beam prior to injection into the wigglers is possible via the 
straight-ahead path (upper path on Figure 4-1). The low power 
beam dumps will permit early tuning of the oscillator 
accelerator, both during APLE construction and afterward when 
in operation. 

The amplified light is produced via the systems two 


undulators. The oscillator wiggler will be designed by STI 
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Optronics, Inc. [32]. It has a period A=2.36 cm and a length 
L=2.36 m, thereby consisting of N=100 periods. The mirror 
Separation Gistanee, Ss, 1s  22eene The APLE oscillator is 
designed to produce a relatively low-power seed laser bean, 
typically 100 watts. This power can be easily found by taking 
the oscillator potential of 17° MV with the currene ores, 
to gét a power of “PeMW “(at e25+ Gury). The extraction 
efficiency of 0.01% yields the laser power of 100 W. This 
laser beam passes through the amplifier concurrently with the 
electron beam to produce the 100 kW laser beam. The seeaq 
laser pulses produced by the oscillator are directed to the 
amplifier wiggler by a series of optics. 

The amplifier wiggler is L=9.96 m long, and consists of 
N=256 periods which are each 4=3.89 cm in length; design 
extraction efficiency is 5%. Both wigglers have a tapering 
capability. The optical beam pulses must properly overlap the 
electron beam pulses to maintain coupling in the amplifier 
section {12}. The required overlap of the short electron and 
optical pulses is achieved by control of the path length and 
direction of the optical beam using adjustable mirror mounts, 
two in the resonator and one in the beam train between the 
oscillator and the amplifier. The laser beam then propagates 
from the amplifier in a 27 m long vacuum pipe to an optical 
diagnostic area [32]. 

Heating of the oscillator mirrors can be a problem; 


therefore, it is important to know the laser beam intensity on 
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the mirror surface. A critical motivation for using the MOPA 
COmmiguration is to circumvent mirrorm heating problems through 
use of the high-gain amplifier, which has no mirrors. The 


first step is to solve for the laser spot size at the 


resonator mirror, w?=w,(1+z7/z/), where w is the spot radius 


aecePpositvon 2,9 is the distance to the mirror from the mode 


waist, and Z,) is the Rayleigh range. The Rayleigh range is 
defined as z=nw,;/A, and is 29=L/2=118 cm for the oscillator 


resonator. Then, the spot area on the mirror is 11.2 cm. 


With a laser power of 100 W, the intensity is 8.9 W/cm? 
deposited on the mirrors. This intensity is quite modest and 
will not cause any heating problems. 

One of the major improvements of APLE over Boeing’s 
previous experiments should be improved beam quality [34]. 
The successful operation of any FEL is critically dependent 
upon the quality of the electron beam. As discussed in the 
previous chapter, beam quality is generally described by two 
quantities: the beam energy spread, Ay/y, and the beam 
emittance, e€. The beam emittance, €, 1s established by the 
injector {35}. The high beam quality must be maintained as 
the beam is accelerated and transported to the wiggler. 
Boeing’s previous FEL [36] exhibited the following: Ay/y=.005 
(FWHM) at 250 Amps, and €,=ye=100-120n mm-mrad. The design 
values for APLE are: Ay/y=.004 (FWHM) at 450 Amps, and ¢€,=40n 


mm-mrad (FWHM). In order to maximize beam quality, Boeing 


Sy, 


chose the photoelectric cathode injector over the thermionic 
oun. APLE will be the first to use a laser-driven, 
photocathode injector in high-duty operation [32]. 

Recall that a few dimensionless variables can summarize 
recurring physical properties in the FEL. Table 4-2 shows 


some of the important FEL dimensionless parameters. 


TABLE 4-2. APLE DIMENSIONLESS PARAMETERS. 


jJ, current 250 197000 
density 

K  pundudiaror O23 dee 
parameter 
number of 
periods 


0.0094 


Pa 
— 







Pot Lee ace ee 
weak fields, 
Single-mode 









Pe Lv lange tractor 0.04 
strong fields, 
multi-mode 


Og, gaussian 20 
distribution 
energy spread 
—e 











Og, exponential 
Gust risue1on 
energy spread 
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V. THE APLE OSCILLATOR 


A. DIFFRACTION EFFECTS 

In this chapter, numerical simulations are used to 
investigate important features of the APLE oscillator. These 
Simulations use the parabolic wave equation [12] to solve for 
the natural diffraction and transverse mode distortion in the 
optical wavefront. Only one field site in the longitudinal 
dimension is followed, but the simulation output includes 
evaluation of the FEL efficiency, gain, and the transverse 
optical mode - electron interaction. 

One of the initial steps in the investigation is to find 
Piet ial Fee phase velocity, v5, fom maximum gain. The 
initial phase velocity determines the FEL gain at the end of 
ene wundulator [12]. A gain spectrum is obtained by 
calculating the gain as a function of phase velocity. The 
peak gain is then selected directly from the gain spectrum, 
and the corresponding value of v, is used in subsequent 
numerical simulations. 

The APLE dimensionless parameters are used as inputs to 
Efe wmmumermrecat SOlution Of the @gain spectrum. The 
dimensionless current density is j=250, and the radius of the 
electron beam is o,=r,(m/L A)'*=0.28. The beam quality is 
accounted for by using the distributions og=1.5 and o,=1.0. 


The dimensionless Rayleigh length, 2)=Z,)/L=0.5, is found by 
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dividing the dimensional Rayleigh length by the undulator 
length LZ. The position of the optical mode waist along the 
undulator is described by t,=0.5 when centered at the middle 
of the undulator at £/2. The initial omensionilescme seca 
field strength is estimated as 4,=0.01 in the weak field 
regime. The electron beam traverses through N=100 periods 
dumang a single tripeeneougn™ = emia velce or: 

Figure 5-1 shows the spectrum of final gains, G(v)), 
plotted against the initial phase velocity v). The peak gain 


of G=23.5 occurs at v,=5, and the spectrum is nearly symmetric 


x** FEL 3d simulation, single-pass gain *** 


3=250 Gg=l.5 © =O .28 Og=1 


a=0.01 N=100 





Figure 5-1. Gain spectrum of the APLE oscillator. 
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about vo=5. If the gain bandwidth is defined by the range of 
phase velocities over which the gain is reduced by =1l/e, then 
AVo=2TT. Wnen Snatural edie nmactlonmrs sume luded, the optical 
phase changes according to $(t)=-t/z,j=-2t along the undulator, 
angedaves a shift in resogance by avj=1/z,=2. [12]. The energy 
spread shifts the gain spectrum to higher values of v, so that 
Elemecameqaim Should occur Close tO )=Se8as is seen. The 
gain spectrum shows how the gain changes when the optical 
wavelength varies through av,=2mNaA/A around the resonant 
wavelength, or when the electron beam energy varies through 
AV>=4mNay/y around the resonant energy. 

The FEL wavefront evolution as t goes from 0 to 1 along 
the N=100 period undulator is illustrated in Figure 5-2 via 
the results of a numerical simulation. The input parameters 
Cee clita to those in the gain spectrum simulation. The 
C=rcCaomOort DOM cChesGy Spredgeio.) and semittance (cg) ane 
ieeoimeedenor sand we use the additional parameter v)=5. 

ihe evolution One ei @peical mode amplitude, 
|a(x,t)|(upper-left), is shown as an intensity/contour plot. 
iMcmic imum tela iseblack andsthe zero field as white; the 
white contour line represents the field at 50% of the maximum 
amplitude shown on the right. The grey scale associated with 
these plots is shown in the upper-right window. The evolution 
shows that the peak field, in the center of the optical mode, 
grows steadily through the undulator. The upper-center frame 


shows the mode, |a(x,y)|, at the end of the undulator; x and 


43 


y axes are normalized by (LA/m)', where L=236 cm and 210.6 


um. The evolution of the bunching current, “o(77 7) sqm cee. 


2\||a(x,T) | 







a ’ 


0.2 
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Figure 5-2. Numerical simulation of the APLE oscillator 


wavefronts in weak fields. 


left), shows the amount of bunching in the electron beam along 
the length of the undulacorm: There is no bunching at 
the beginning of theP undulator ny (Tr=0)— buem@iunere 25 ecleaa 
bunching at the’end of the undulaton (ral). ~The midale-cencer 
frame shows the electron bunching at the end of the undulator 
o(x,V). Recall that themtusling Gaeeemas: ols schesmatiG jot 
the electron beam area, =m(o,)* where o,=0.28, to that of the 
optical mode. The focussing of the optical mode increases the 
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value of F, and serves to increase the gain. The electron 
Phase. velocity .distrisweren, fivemeilower-left), shows 
negligible distortion of the electron beam in the weak optical 
field. The final electron phase-space distribution, (¢€,v) 
(lower-center), shows a small amount of electron bunching at 
the end of the undulator. The spread in phase velocities due 
eomene wet ialespread o.=1 and og=I5 1s clearly visible. The 
plot of the optical field phase evolution at the center of the 
mode, 9$(0,t) (on the right),. indicates an almost linear 
decrease in phase. This is caused by natural diffraction 
Soeiemcwonase Shitt Of @(ej=-t/z,=-2e over the length of the 
unG@ulator . The power evolution, P(t), shows a progressive 
increase in power as bunching develops. The gain, G(t), also 
increases with the increased bunching of electrons; the final 
Gain@es C.=25., 

Figure 5-3 gives the result of a simulation with nearly no 
current (j=0.001) in order to show the laser mode, |a(x,y)|, 
at the end of the undulator due to natural diffraction. The 


Vavet Eat 15 sinitaalized= to have "Rayleigh length zas0.5 


Preonedminimene mluaqle of thesundulator at t.-0.5, The mode 
walst radius 1s Waatdz,))  ae0. 7 07. At the beginning and end of 
the imGuhalor , the mode radius is given by 


w=w,(1+0.52/z¢)1/2=/2w,=1 as is seen in the figure. The mode 


in figure 5-2 is clearly smaller with a radius of w=0.4, as 
the high current electron beam has focussed the light inward 
toward the beam. 
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0 6 x 


Figure 5-3. Optical mode at the end of the undulator, when 
only governed by diffraction. 


Figure 5-4 is another numerical simulation of the FEL 


wavefronts, but with a different input parameter. The 





0 15 1-nx/2 ¢ 3n/20 ie a 


Figure 5-4. APLE oscillator wavefronts with an exponential 
energy spread, and a Z,=0.28. 
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Rayleigh length z,=(1/12)'/*=0.28 vice z)=0.5. The smaller 
Rayleigh length of z,=(1/12)' focusses the optical wavefront 
around the electron more tightly all along the undulator, and 
increases coupling. The overall result is a significantly 
higher final gain of G,=28. 

We now investigate the oscillator wavefronts through 
successive passes of the undulator. This numerical simulation 
includes inputs similar to those in Figure 5-2, with some 
added parameters. The FEL system simulated now includes 
mirrors. The dimensionless mirror radius r,=2.4 is normalized 
by mi a). The dimensionless radius of curvature of the 
Niger ieees «normalized by L. Figure 5-5 shows the 
results of this simulation for n=32 passes through the 
undulator. The mirror surface near the beginning of the 
undulator is taken to be perfectly reflecting. The mirror 
surface near the end of the undulator is given a loss 
described by Q=5, defined so that over many passes n, the 


-n/Q Each mirror radius is 


power loss would be given by «e 
determined so that 1% of the optical power in the fundamental 
fecdewgoecema Ound the Mirror, as indicated by edge loss. The 
radius of curvature of each mirror is determined by the 
Rayleigh length of the optical mode, z,=0.5, and the position 
See oemncde waist, CC —0.5. SThe@distance between the mirrors 
is taken to be twice the undulator length with the undulator 


eéneered between the mirrors. This resonator length is much 


less than anticipated in the experiment, but allows a 
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numerical solution and represents the appropriate Rayleigh 


length. 


5{/a(x,n) | 40) |a(x,y) | 
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0) n S272 26 3n/20 n 32 
Figure 5-5. APLE oscillator wavefront evolution over n=100 
passes. 


The evolution of the optical mode amplitude (upper-left), 
shows that the field amplitude grows with each successive pass 
down the undulator until reaching a maximum amplitude of 
ja(x,n)|=40. The upper-center frame is an “end-on" view of 
the optical mode at the end of the undulator on the final 
pass. The middle window depicts the optical field amplitude, 
la(x,t)|, along the length of the resonator from mirror to 
mirror on the final pass. Notice that the largest field 
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amplitude is near the center of the undulator at the focus of 
the mode waist near t,=0.5. The lower-left window shows the 
evolution of the electron phase velocity distribution, f(v,n), 
over successive passes. The electron phase-space evolution 
(¢,v) at the end of the final pass (lower-center), shows 
Significant electron bunching in closed phase-space paths and 
strong optical fields. The power evolution with each pass, 
P(n), shows a steady increase in power followed by relatively 
Steady power at saturation. This is also reflected in the 


plot vemeagaim™ lower-right), Gin). 


B. PULSE EFFECTS 

In this section, we use numerical simulations to 
investigate the multimode pulse effects of the APLE 
oscillator. We begin by observing the electron phase-space 
evolution, and then investigate the FEL pulse evolution. 

Figure 5-6 shows the result of solving the pendulum and 
wave equations numerically with 2000 sample electrons in a 
strong optical field a,=4n*=40 and for low current of jF=10 
where j=250 and estimating F=0.04. The electrons are started 
Vearieecmecandom Gaussian spread ‘corresponding to o,=1 and 
emittance corresponding to og=1.5 about the phase velocity 
Vo=5- The final phase-space positions of the electrons are 
drawn in (¢,v). The final positions of the electrons show 
overbunching, with many electrons trapped in the closed orbits 


of phase-space. Some of the electrons show an increase in 
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energy, or phase velocity v. These electrons absorbed energy 


from the optical mode, resulting in a decrease in gain at the 


**xk* FEL Phase Space Evolution *** 
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Figure 5-6. APLE oscillator strong-field phase-space 
evolution with low current. 
end of the undulator. This decrease is evident in the 


corresponding curve to the right. 

Figure 5-7 shows the APLE oscillator pulse evolution over 
n=800 passes with no initial field (a =0), starting at 
resonance v,=0. The window width w=30 corresponds to thirty 
slippage distances along the beam, z is the longitudinal axis. 
The dimensionless current density is jF=10. The dimensional 
electron pulse length is normalized by NA to get the 
dimensionless electron "pulse) Jengeno ase The 
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desynchronism (see section III. G.), defined by d=0.0001, 


indicates a very small change in the distance between mirrors 


tn FREePalse hyotutcion *~**x-* 
Oo =17 d=0.0001 Og=1 a5 
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Figure 5-7. APLE oscillator pulse evolution. 


of aS=0.1 um. A small random phase, 6€=10%, is added to each 
of the electron phases to represent shot noise. The shot 
noise provides spontaneous emission for starting the FEL 
Goch our ine electron beam energy spread is o2=1.5>1 to 
account for both exponential and Gaussian distributions. 
However, when the oscillator saturates and the fields become 
strong, the shape of the distribution is of little importance. 


The factor Q=5 determines the resonator loss on each pass 
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through the undulator. The undulator consists of N=100 
pebiedear The evolution of the optical field amplitude 
|a(z,n)| (middle-left), shows unstable oscillation ana Ee 
trapped-particle instability. The optical pulse length is =10 
(FWHM) slippage distances. The window above shows the final 
optical field, where the modulation of the field envelope is 
clearly evident. The modulation is caused by the trapped- 
particle instability. The grey scale above |a(z,n)| shows the 
peak field amplitude, |a(z,n)|=45, in white, and zero field in 
black with one contour. Recall that a,=40n is required for 
one oscillation of trapped-particles. The power P(n) 
increases as sidebands appear in the power spectrum, P(v,n) 
(middle-center). The final power spectrum, P(v) (top-center), 
has a significant sideband with optical power comparable to 
the fundamental. Again, the trapped-particle instability is 
evident. The pointed tick-mark at the top of the final power 
spectrum indicates the central wavelength of the initial 
radiation at resonance [12]. The rectangular tick-mark 
indicates the center of the final power spectrum. For 
reference, the weak=-field gain SpeCerun, sG(y), tombe lovee 
power spectrum on the same scale. It is calculated separately 
for JjF=10 with ne pulses 117). The evolution of the 
electron’s phase-velocity distribution after each pass, f(v,n) 
(middle-right), indicates a relatively wide electron 
distribution with a narrow base. The final phase-velocity 


spectrum of the electrons in the beam, f(v) is shown at the 


S)Z 


EOD-—cpight. The pointed tick-mark indicates the electron 
beam’s initial phase-velocity at resonance. The position of 
Peewee lectron Pulse, 7(z-t)F (lower-left), is indicated in dark 
iaey dteecere beginning of the undulator, c=0, and is light grey 
at the end of the undulator, t=1. The unbunched electron 
pulse starts at t=0 and slips back to its position at t=1 on 
each pass, as bunching develops. The lower-right window shows 
the optical power, P(n), after each pass. The optical power 
appears to reach steady state after n=300 passes. These 
observations correspond to the evolution of the optical field 
amplitude in the middle-left window. 

Figure 5-8 shows the same oscillator as in Figure 5-7, 
but with a desynchronism of d=0.01. The electron spectrum 
1s broader, this is especially evident after =n=250 passes 
whem me@npeselectronsegive up energy to the optical field. The 
optical pulse is longer, and forms coincidentally with the 
initial spread in the electron spectrum at =n=50 passes. The 
lonG@emepepticaimepulse is beneficial, because more of the 
electron pulse (lower-left) can interact with the optical 
field as the pulse slowly slips back from t=0 to t=1l. This 
serves to increase the gain. The optical field evolution and 
mae ssspOwer €vOlution still show the trapped-particle 
instability. The oscillator exhibits relatively stable 
operation. 

Figure 5-9 shows the oscillator with a slightly larger 


desynchronism, d=0.03. Here, the electron spectrum is 
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xxkxk FEL Pulse Evolution **** 
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Figure 5-8. APLE oscillator with d=0.01. 
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Figure 5-9. APLE oscillator with d=0.03. 


a5 


broadest, extending furthepmem the -v direction ne soon 
that more electrons are giving up energy to the optical field. 
The optical pulse is also longer, forming coincident with the 
onset of the broad electron spectrum at =n=100 passes. The 
trapped-particle instability still exists. The power may be 
unstable because of chaotic optical modes. 

Figure 5-10 shows the same oscillator with a still larger 
desynchronism, d=0.08. The broad electron distribution still 
exists, although it is slightly narrower than that of the 
previous figures. The optical pulse length is about the same 


as in Figure 5-9, but its amplitude has decreased. Notice 
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Figure 5-10. APLE oscillator with d=0.08. 
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that the power spectrum has a smaller sideband. The optical 
power is also less than that exhibited with d=0.03. However, 
the FEL operation appears to be stable. 

Figure 5-11 shows the same oscillator as in Figure 5-10, 
but with a desynchronism d=0.1. The trapped-particle 
instability has decreased, as evidenced by the optical field 
evolution and the power evolution. The electron spectrum is 
Slightly narrower than shown in Figure 5-10, and the optical 
pulse is shorter. These factors combine to cause a decrease 
in power. 

We cannot explore larger values of desynchronism because 
of n@mertca® limitations ® but the trends indicated in the 
previous examples continue. As desynchronism increases 
further, the electron beam distribution is made smaller and 
the optical pulse is) longer. POSE Oumechemelieces are 


desirable for the amplifier performance. 
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Figure 5-11. APLE oscillator with d=0.1. 
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Vi. THE APLE AMPLIFIER 


A. DIFFRACTION EFFECTS 

In this chapter, we use numerical simulations to 
investigate important features of the APLE amplifier. The 
method used to evaluate diffraction effects is similar to that 
used in the previous chapter. However, the pulse simulations 
evaluate a single-pass (t=0 71) vice multiple passes, because 
the light and electrons traverse the amplifier only once. The 
evaluation of the amplifier is also different because it 


includes taper. With tapering included, the pendulum equation 
can be written as t = §6 ( |a|-a.) +jalcos (+o) where optical field 


amplitude at saturation = 2 WI a 6 is a step function, and 


6 represents taper [12]. 

The first step in this investigation is to follow the FEL 
wavefront evolution as t goes from 0 to 1 along the N=256 
period undulator, as illustrated in Figure 6-1. iive 
dimensionless current density is j=197000, and the radius of 
the electron beam is o,=0.09. Rhiese Pie ial optical. fie) 
amplitude a4,=18. A gain spectrum simulation (Similar to 
Figucve 5-1} showed that peak gain is at v,)=15. The wavefront 
Sl amonerS GDuMeat Vo=l5 tor peak gain in weak fields prior 


2/5 


to Saturation at field value 4a,=2(j/2) The beam quality 


is accounted for by using the distributions oOg=7 and o,=6. 
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The dimensionless Rayleigh length 2,=0.8. The optical mode is 
focussed at the beginning of the undulator at t,=0. Moke 


amount of taper is indicated by 6=-4nNK* aK/K(1+K*)=400n. 


4)|a(x,T) | lja(x,y) | 
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Figure 6-1. APLE amplifier wavefronts in strong fields. 


The evolution of the optical mode amplitude, |a(x,T) | 
(upper-left), shows that the peak field grows steadily through 
the undulator. Without focussing by the intense electron 


beam, the size of the optical mode at the end of the undulator 


would be w=/Z)(1+1/z5)1/2=1.43 where 29=0.8. Instead, the 
optical mode radius 1s reduced to w=0.3. We obtain the 


filling factor by taking the ratiomomeeresc (|CCeron Dealeasea, 
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Gx;yy (middle-center), towthe optical mode area, |a(x,y)| 
(upper-center). [his @patio™resu lism eas filling. factor 
F=0.10. The focussing of the optical mode increases the value 
of F, and serves to increase the gain. The evolution of the 
bunching current, o(x,t) (middle-left), shows there is no 
pomening at the beginning Of she wadulator #r=0) ,ebut@there is 
clear bunching at the end of the undulator (t=1). The 
electron phase-velocity distribution, f(v,t) (lower-left), 
shows that about half the electrons become untrapped after 
t>0.25 when the optical field reaches the saturation level 
a.=2(j/2)*. Experience indicates that when in strong optical 
fields, the tapered FEL operates best when *=50% of the 
electrons remain trapped. Increasing the taper serves to 
untrap more electrons. The final electron phase-space 
distribution, (¢€,v) (lower-center), shows a large amount of 
electron bunching at the end of the undulator. Note that =50% 
of the electrons are trapped. The spread in phase velocities 
due Beemenece wiatral Spreadgo.~6 “and o,~/ is clearly visible. 
The plot of the optical field phase evolution at the center of 
the@meee ww er0, tc) (On mtne right), indicates anyincrease in 
phase followed by a steady decrease. The power evolution, 
P(t), shows a progressive increase in power as bunching 
develops. The gain, G(t), also increases with the increased 
Durie mgmemeelectrons;ethe finaleagqain is 9¢.=25000. 

It can be shown that =10% of the optical power, an average 


of =~10 kW, strikes the inside surface of the beam pipe. Over 
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heating of the beam pipe may occur; therefore, this area 


should be investigated further. 


B. PULSE EFFECTS 

Figure 6-2 shows the optimum APLE amplifier pulse 
evolution over one pass with an initial field a,=18, starting 
at v,)=15. The window width w=4 corresponds to four slippage 
distances along the beam, z is the longitudinal axis. The 
dimensionless current density is jF=20000, assuming F=0.10. 
The dimensionless electwon pulse, lengths )o-=2) (iii uae 


undulator taper is 5=-4nNK°sK/K(1+K*)=250n, corresponding to 
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Figure 6=2. APLE ampiet Ben pulse sey Gee on wap aD. 
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aetaper of about 40% mm the undulator field strength. The 
optical pulse width, o,=6 (FWHM), is found by taking the 
oscillator optical pulse length and normalizing it for the 
amplifier. The undulator consists of N=256 periods. The 
evolution of the electron’s phase-velocity distribution, 
f(v,t) (middle-right), must be viewed differently than in the 
case of the oscillator, since the amplifier is tapered (see 
Chapter III. F.). Recall that tapering maintains the optical 
field resonant with the electrons as they lose energy. 
Electrons trapped near resonance v=0 actually lose energy as 
they proceed through the undulator. The distribution shows 
electrons becoming untrapped about 40% of the way down the 
undulator when |a|>a,. The trapped electrons, maintaining a 
phase velocity v=0 down the undulator, lose energy to the 
optical field thereby amplifying the light. The final phase- 
velocity spectrum of the electrons in the beam, f(v), 
indicates that =50% of the electrons are trapped. The 
evolution of the optical field amplitude |a(z,t)| (middle- 
left), shows growth coincident with the untrapping of 
electrons. The grey scale above |a(z,n)| shows the peak field 
amplitude, |a(z,t)|=4436. The power spectrum, P(v,T) 
(middle-center), is narrow. The final power spectrum, P(v), 
has no sidebands; the trapped-particle instability is not 
present. Recall that the position of the electron pulse, 
j(z-t)F (lower-left), 1S indicated in dark grey at the 


beginning of the undulator, t=0, and is light grey at the end 
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Of the® Unee@tavon =e c=ie It is evident that there is 
Significant optical field - electron pulse interaction.» The 
lower-right window shows the optical power evolution through 
the undulater, = +(e). The optical power appears to grow 
steadily after electrons become untrapped at 1t=0.4. The 
efficiency 7=8.8%. 

The use of a short optical pulse is an option being 
considered by researchers at Boeing. Therefore, we now 
investigate the amplifier’s performance with a shorter optical 
pulse, now the same length as the electron pulse. Figure 6-3 
shows the same amplifier as in Figure 6-2, but with an optical 
pulse o,=2. The electron spectrum indicates slight under 
trapping of electrons. Notice that the final power has 


decreased. The efficiency is also less, n=7.6%. 
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Figure 6-3. APLE amplifier pulse evoltion with o,=2. 
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Vil. CONCLUSIONS 


The demonstration of high average power is one of the 
important directions that FEL research is heading. When 
coupled with the FEL’s demonstrated tunability and 
reliability, high average power will create a highly desirable 
directed energy system. The FEL’s quick response time and 
"infinite magazine" also make it a candidate for military 
applications. Military applications include shipboard, land- 
based and airborne FELs. The advantages of directed energy 
and the FEL weapon specifically, were discussed herein. The 
FEL is likely the weapon of choice for the next century, and 
could have significant impact in many warfare areas, including 
Theater Ballistic Missile Defense (TBMD). The above features 
make the FEL @attragtave for sei1emeiriic ana Inedstria 
applications as well. 

The proposed Boeing Average Power Laser Experiment (APLE) 
will demonstrate the highest FEL average power to date (100 
kW); completion is scheduled for 1996. The APLE design 
appears feasible, and can result in successful attainment of 
the design goal. 

The research addressed in Chapters V and VI provide some 
insight into the experiment. A relatively large desynchronism 
is beneficial to the oscillator. This results in stable 


performance and a long optical pulse. The long pulse, when 
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compared with a short pulse, optimizes the performance of the 
amplifier, resulting in both higher efficiency and higher 


power. 
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